Abstract. We include the tidal instability due to the 
3:1 resonance in the disc instability model developed by 
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Hameury et al. (1998) and modified by Buat-Menard et 
al. (2001). We confirm earlier results by Osaki (1989) that 
the tidal instability can account for the SU UMa light 
curves. We show that in ultra-low mass ratio systems such 
as ER UMa stars and WZ Sge stars, the superout burst 
ends while the disc is still eccentric, as proposed by Hel- 
lier (2001). However, since the disc shrinks rapidly once 



a cooling wave has started, the eccentricity should stop 
shortly after the end of a superoutburst. This result dis- 
agrees with the suggestion by Hellier that decoupling the 
thermal and tidal instability in the TTI model can account 
for late superhumps and echo outbursts in ultra-low mass 
ratio systems. We propose instead that ER UMa short 
supercycles can be explained either by the alternation of 
narrow and wide outbursts similar to those occurring in 
SSCyg, or by the effects of irradiation (Hameury et al. 
200C). In both cases, we predict that superhumps should 
be permanent, which is suggested by observations (Gao 
et al. |1999P . We can also reproduce light curves similar to 
those of EG Cnc, varying the mass transfer rate in a TTI 
model including both irradiation and the presence of an 
inner hole in the disc. 
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1. Introduction 

Dwarf novae are cataclysmic variables (CVs) which ex- 



hibit 2-6 mag recurrent outbursts (see Warner |1995 for an 
encyclopaedic review). It is generally accepted that these 
outbursts are due to a thermal- viscous instability that oc- 
curs in an accretion disc in which the viscosity is given by 
the a-prescription (Shakura & Sunyaev 1973Q . The insta- 
bility arises when hydrogen becomes partially ionized and 
the opacities vary steeply with temperature; this occurs 
when the disc material reaches temperatures of order of 



8000 K (see Lasota 2001 for a review of the model). A num- 
ber of characteristics of many light curves (characteristic 
time scales, amplitudes, ...) are reproduced by the disc in- 
stability model (DIM) provided the viscosity parameter a 
is different in quiescence and in outburst. However several 
subclasses of dwarf novae show specific features that are 
not easy to explain with the standard version of the DIM. 

SU UMa stars are defined as dwarf novae in which nor- 
mal outbursts cycles are regularly interrupted by anoma- 
lously bright and long outbursts called superoutbursts. 
These outbursts are 0.7 magnitudes brighter and last 
5 to 10 times longer than the normal ones. The mean 
time T s between two consecutive superoutbursts defines 
the period of a supercycle T s ; in most SU UMa stars, 
130 < T s < 350 days. There are however two extreme 
exceptions: ER UMa stars have anomalously short su- 
percycles (between 19 and 45 days), while WZ Sge stars 
show only superoutbursts separated by very long quies- 
cence periods that can reach 30 years. In addition, the 
light curve is modulated during superoutbursts at a period 
slightly longer than the orbital period. These so-called su- 
perhumps usually disappear short ly after the decline of the 
superout burst. Whitehurst (198S) explained these super- 
humps by the precession of a distorted disc. The name 'su- 
peroutburst' is reserved for dwarf nova eruptions in which 
a super hump is present. This property differentiates su- 
peroutbursts from wide outbursts observed in many dwarf 
novae; van Paradijs ( 1983 ) suggested that superoutbursts 
are just wider outbursts in the bimodal outburst-length 
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distribution. In addition, a 1985 outburst of U Gem had 
all the properties of a superoutburst (of WZ Sge stars) 
except for the superhump (see Kuulkers et al. 199£ and 
references therein). This could suggest that the superout- 
burst phenomenon is independent of the disc distortion. 
We discuss this possibility in the Section 4. 

As almost all SU UMa stars have periods below the 
period gap of CVs, they all have low secondary to pri- 
mary mass ratio q, and the primary Roche lobe is large. 
The disc can therefore extend to relatively large radii, and 
this led Osaki ( |1989| ) to propose the thermal-tidal insta- 
bility (TTI), in which superoutbursts are caused by the 
disc reaching the 3:1 resonance radius, as a model for SU 
UMa stars. SPH simulations by Whitehurst (1988) showed 
that the disc becomes eccentric and precesses when its ra- 
dius reaches the 3:1 resonance radius, but it must be kept 
in mind that the hydrodynamical simulations of Stehle 
( 1999| ) and Kornet & Rozyczka ( 2000| ) did not show such 
an effect. 

Detailed observations of SU UMa stars point out sev- 
eral limitations of the TTI model. The mass transfer rate 
from the secondary is observed to be enhanced during 
outburst; this increase, possibly due to irradiation, is not 
taken into account, although it can drastically alter the 
light curves (Hameury et al. [2000 ). 

The TTI model cannot easily reproduce the properties 
of ER UMa and WZ Sge subclasses. In order to obtain 
the very short ER UMa supercycles (19 to 45 days), one 
must assume that the tidal instability shuts off prema- 
turely (Osaki 1995a ) for very small mass ratios. The very 
long quiescence time of WZ Sge systems would on the 
other hand require a very small viscosity in quiescence 
(Smak |l993| ). 

Occasionally, some systems show echo outbursts at 
the end of a superoutburst: several consecutive small out- 
bursts with very short recurrence time are triggered be- 
fore the disc returns to quiescence. During these echo- 
outbursts, superhumps are still present. These so-called 
late superhumps (LS) have also been observed in ER UMa 
stars. However, the TTI model simulations predict that 
the end of a superoutburst should coincide with the ces- 
sation of the tidal instability. 
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These different problems have lead several authors to 
investigate alternative possibilities. Hameury ct al. (2000) 
have shown that if one takes into account the disc and 
secondary illumination in the standard DIM, one is able to 
reproduce light curves reminiscent of those of systems such 



as RZ LMi or EG Cnc. Recently Hellier fl2001| ) suggested 
that some difficulties could be solved in the framework 
of a slightly modified TTI model; he proposed that the 
superoutburst must end before the end of eccentricity in 
low mass ratio systems, thus giving an explanation for LS. 
We have included in our DIM model (Hameury et al. 
Buat-Menard et al. 2001, hereafter Paper I) an en- 



1998 



hanced tidal torque when the disc radius exceeds the 3:1 
resonance radius. We do reproduce the earlier results of 



Osaki ( 1989 ), obtained using a simplified model. We show 
that in a number of cooling wave can end a super- 

outburst while the disc is still eccentric, thereby confirm- 
ing the proposal of Hellier ( 2001 ). However, this by itself 
is not sufficient to reproduce the properties of ER UMa or 
WZ Sge systems, and additional effects such as illumina- 
tion and mass transfer variations must be included. 
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Fig. 1. Superoutburst cycle in the TTI model for SU UMa 
with Mtr = 10 16 g.s~ 1 . Top panel: disc radius; intermediate 
panel: disc mass; bottom panel: disc visual magnitude 



2. The TTI model 

The equation of angular momentum conservation in a disc 



is (Hameury et al. 1998) 
<9£ Id, 
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where S is the surface column density, M tr is the rate at 
which mass is incorporated into the disc at radius r, v r 
is the radial velocity in the disc, j — (GMir) 1 / 2 is the 
specific angular momentum of material at radius r in the 
disc, £1k = (GMi/r 3 ) 1 / 2 is the Keplerian angular velocity 
(where Mi is the primary mass), v is the kinematic vis- 
cosity coefficient, and ji the specific angular momentum 
of the material transferred from the secondary. T t id is the 
torque due to the tidal forces, for which we use the for- 



mula of Smak (1984), derived from the determination of 
tidal torques by Papaloizou & Pringle (1977): 

Ttid = ajrvT, (2) 
where oj is the angular velocity of the binary orbital mo- 
tion, a the binary orbital separation and c a numerical 
constant taken so as to give an average disc radius equal 
to some chosen value. 

In the TTI model, c is no longer a constant, but de- 
pends on whether the disc is eccentric or not. Ichikawa et 
al. ( 1993] ) reproduced supercycles assuming that the disc 
becomes eccentric when its radius reaches r^i ~ 0.46a; 
c is then increased by a large factor. Afterward the disc 
starts shrinking and can no longer maintain its eccentric- 
ity when its radius becomes smaller than r cr ito ~ 0.35a, 
at which point c returns to its normal value. We choose 
to take a value Co when the disc is axisymmetric and a 



larger value c\ when the disc is eccentric. This is our Pre- 
scription / for c, quite similar to that of Ichikawa et al. 



(1993) who extended the work of Osaki (1989) by taking 
into account a non zero tidal torque during quiescence. 

We have included this variable tidal torque in the disc 
instability model used in paper I, which is a modified ver- 
sion of the DIM described in Hameury et al. ( [1998 ) that 
includes heating by the stream impact and the tidal torque 
(We have corrected Eq. (1) of paper I, in which the heating 
rate was overestimated by a factor ~ 2). In the following, 
we also assume that a = a co id = 0.04 in quiescence and 
ct = c^hot = 0.2 in outburst. The inner radius is fixed at 
r in = 10 9 cm. 

We obtain the cycle shown in Fig. [|. After a super- 
outburst, the disc mass and radius are small. The disc 
mass slowly increases, and several normal outbursts are 
triggered due to the thermal instability. During each out- 
burst, only a small amount of mass is accreted by the 
white dwarf, so that the average disc mass increases; a 
significant disc expansion occurs during the rise to max- 
imum, when the heat front reaches the outer edge of the 
disc. The disk then shrinks when a cooling front brings the 
disc back to a cool state, causing a significant reduction of 
the outward angular momentum flow. After the end of an 
outburst, the disc radius very slowly decreases as a result 
of mass addition and viscous diffusion, but the net balance 
is an increase of the radius from one outburst to the next 
one. Eventually, r out becomes larger than r^-x during an 
outburst and the tidal torque is increased, as well as the 
associated energy dissipation. The disc is then maintained 
in the hot state until most of the mass has been accreted. 
During this phase the radius decreases until r < r cr i t o- 
The disc returns to a circular shape, the anomalous tidal 
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dissipation stops, and a cooling wave starts from the disc 
outer edge, bringing the whole disc in a cool state. 

2.1. Radial dependence of c 

There is, however, no compelling reason to assume that 
the whole disc is affected by the development of the tidal 
instability. We therefore use another prescription for c 
(hereafter Prescription II). We assume that when the 
disc is circular, c = Co; when eccentric, c(r) = cj for 
r > 0.9 r cr ito ancl c = Co for r < 0.8 r cr ito- A linear 
interpolation is assumed between 0.8 r cr ;to and 0.9 r cr ito- 
We compare and discuss the results obtained using both 



Table 1. Parameters used for SU UMa subclasses simu- 
lations. Mi is the secondary mass and r C j rc the circulation 
radius. 



prescriptions in 3.1 



2.2. Time delays 

Temporal variations of c(r) are not instantaneous. Obser- 
vations show that normal superhumps appear within a few 
of days after the onset of a superout burst, and that they 
develop within a day, so that it does not take long for 
the disc to become eccentric (Semeniuk 1980 ). However, 
it takes the disc at least a dynamical time to change its 
geometrical shape; we assumed that, when the tidal insta- 
bility sets in, c increases linearly on a timescale t r = 2 10 4 
s. Note that in WZ Sge systems, the early superhumps 
that first appear seem to be of a different nature from the 



normal superhumps; Osaki & Meyer (2002) suggested that 
they could be related to the 2:1 resonance. 

Similarly, the disc must dissipate its excess energy in 
order to return to a circular shape. The relative energy dif- 
ference between a circular and an elliptic orbit with eccen- 
tricity e with the same angular momentum is AE/ E c - lrc = 
e 2 . As the dissipation at radius r is D ~ GM'Ev/r 3 , the 
time needed for the disc to return to a normal state is: 



tr. 



AE 



e 2 U 



(3) 



D 2v 

where t v is the viscous time in the outer part of the disc. 
When the disc reverts to its normal state, the cooling wave 
has already started or will start with the end of tidal insta- 
bility (see section ||). t v therefore refers to the cool state, 
and is of the order of the recurrence time, since the out- 
bursts obtained here are of the inside-out type. Typical 
values are tens of days. Numerical simulations of the de- 
velopment of the tidal instability show that e is of order of 
0.1-0.3 (Murray 1996|) , which is confirmed by observations 
(Patterson et al. 200C), and £ c _ c should be a few percent 
of t v . We then assume that c linearly decreases when the 
disc returns to a circular shape on a time if = 2 10 5 s. 

We use three different sets of parameters that are rep- 
resentative for SU UMa, a typical ER UMa star and a 
typical WZ Sge star respectively, shown in Table |l|. r c ; rc is 
the circulation radius at which a particle leaving the La- 
grangian point would stay in circular orbit if there were 
no accretion disc (see Paper I for details). One should 
note that r c ; rc is calculated taking into account the grav- 





SU UMa 


ULMR Stars 






ER UMa type 


WZ Sge type 


Mi/M Q 


0.8 


1.0 


0.6 


M 2 /M Q 


0.15 


0.1 


0.1 


Porb (hr) 


1.83 


1.522 


1.3752 


a/10 10 cm 


5.2 


4.82 


3.876 


r-circ/10 10 cm 


0.91 


1.1 


0.74 



itational potential of the secondary, and therefore differs 
slightly from the Keplerian value. In the following section, 
we will use the term ultra-low mass-ratio (ULMR) stars 
when referring to simulations with ER UMa parameters. 

3. Decoupling the thermal and tidal instabilities 

3.1. A property of ultra-low mass ratio systems ? 

Hellicr ( |200l| ) suggested that the ultra-low mass ratio of 
ER UMa and WZ Sge stars is responsible for their pecu- 
liarities. In their case, the Roche lobe extends far beyond 
the 3:1 resonance radius, and the disc could in principle be 
larger than r^-i- Hellier proposed that regions outside the 
resonance radius suffer a lower tidal dissipation than the 
portions in the 3:1 resonance zone. Then, during a super- 
outburst, heating by the tidal torque will not be as efficient 
as in other SU UMa stars, and a cooling wave could start 
while the disc is still eccentric, and will remain so several 
days after the end of the superout burst. 

This would then account for late superhumps. In ad- 
dition, Hellicr suggested that this model could also ex- 
plain the echo outbursts in WZ Sge stars and the short 
ER UMa supercycle: at the end of a superout burst, the 
disc would be still eccentric, and the tidal dissipation suf- 
ficiently large to trigger short outbursts during which the 
heating wave does not propagate to the outer edge (one 
should add that such outbursts are not only short but also 
have amplitudes lower than outbursts in which the heat- 
ing front reaches the outer disc's edge). ER UMa stars 
would then be systems always in an eccentric state, while 
WZ Sge stars would finally return to complete quiescence 
with circular discs. The difference between the supercycle 
duration of ER UMa and WZ Sge stars would be due to 
higher mass-transfer rates in ER UMa's. 

Fig [j] shows the simulated light curve for the param- 
eters of SU UMa, assuming a factor 20 between c\ and 
Co- The superoutburst ends when the disc radius becomes 
lower than r cr i t o, i.e. the end of the tidal instability causes 
the end of the superoutburst. Fig || shows the simulated 
light curve for the parameters of ER UMa, with the same 
Ci/cq. The mass transfer rate of 10 16 g s" 1 chosen here is 
close to the critical value above which the system is stable 
(within a factor 1.3). In this case, as predicted by Hel- 
lier, the superoutburst ends because the thermal instabil- 
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Fig. 2. Same as fig 1 for ER UMa parameters. 



ity stops, while the tidal instability is still effective. How- 
ever, contrary to Hellier's hypothesis (2001), this phase 



is very short because of the rapid shrinking of the disc. 
Therefore, the duration of the supercycle is much larger 
than observed (there is practically no quiescence phase in 
real ER UMa stars): short outbursts are due to the usual 
thermal-viscous instability. As a consequence, contrary to 
Hellier's assumption, the heating front in short outbursts 
always reaches the outer disc's edge. Mass transfer rates 
closer to the stability limit do not lead to shorter supercy- 
cles, but rather to longer superoutbursts; this is actually 
the reason for which it is impossible to obtain very short 
supercycles in the standard TTI mo del with out changing 
the tidal instability condition (Osaki 1995a). 

We have investigated the influence of the ratio Cx/cq 
and of the Roche size on our results. Fig || compares the 
superoutburst behavior for ci/cq factors of 20, 50 and 80 
for ULMR stars with prescription I. As can be seen, the 
width of the superoutburst strongly depends on c\. Very 
small and very large values of c\ both give short superout- 
burst, either because the tidal heating is too small (small 
ci ) , or because the disc contraction under the effect of the 
tidal torque is too large (large c\ ) . Typical c\ / cq should be 
in the range ^ 20 - 50 in order to reproduce the observed 
superoutburst durations. 

The decoupling of the tidal and thermal instability 
was found for values of ci/cq lower than 80, for both pre- 
scription I and II for ER UMa parameters, thereby con- 
firming the first part of Hellier ( 2001 ) proposal. On the 
other hand, for the parameters of SU UMa, the instability 
stopped because the disc shrinks below 0.35a, except for 
ci/co < 10; this limit is increased to ~ 20 when one uses 
prescription II. 

These results strongly support the idea that the decou- 
pling of thermal and tidal instability is a property of the 
ultra low mass ratio systems. One should however keep in 
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Fig. 3. Radius behavior during superoutburst for ultra- 
low mass ratio stars with M tr = 10 g.B -1 . The dot- 
dashed lines represents the 0.46a (upper) and 0.35a 
(lower) radius. For ci/co factors of 20 (dashed line) and 
50 (solid line), the superoutburst lasts ~ 21 days (~ 17 
days for Ci/cq = 50) and is terminated by a cooling wave. 
For a factor 80 (solid line) it lasts ~ 9 days and ends by 
the disc becoming circular. 



mind that this conclusion depends somewhat on the as- 
sumed radial dependence of c, and on the assumed C\/cq; 
extreme values (less than 10 or more than 50) cause too 
short superoutburst to be triggered but ER UMa stars 
sometimes have short superoutbursts. 

3.2. The superoutburst cooling wave 

As shown in the previous section, decoupling of tidal and 
thermal instabilities is expected in ultra low mass ratio SU 
UMa stars. However, this is not sufficient to explain the 
late superhumps or echo outbursts observed in E R UMa 
and WZ Sge systems, as suggested by Hellier ( 2001 ). When 
a cooling front starts to propagate, the radius decreases 
rapidly (in less than a few days) to its minimum value in 
the supercycle. If this minimum is less than r cr i t o, then 
the disc will return to a circular shape in a time scale 
tf, and one does not expect to observe late superhumps. 
If on the other hand the minimum radius is larger than 
r cr ito, the disc always remains eccentric, the superhumps 
are permanent, but one then loses the very mechanism 
that was supposed to trigger superoutbursts. 

The disc shrinkage at the end of an outburst or a su- 
peroutburst is therefore a key question. The radius evolves 
under the influence of three main effects, as shown in Eq. 
([!]) : (i) the addition of mass with low specific angular mo- 
mentum at the outer edge of the disc and (ii) the tidal 
torque that both act to contract the disc, and (iii) the vis- 
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Fig. 4. Responses of the outer radius of a steady accre- 
tion disc in the hot state to variation by a factor 2 of M tr 
(solid line), c (dashed line) and a (dot-dashed line). The 
parameters vary linearly in 10 4 s. Note that the perturba- 
tion starts at t ~ 0.3 d. 

Table 2. Response of the outer disc radius rd to a de- 
crease of |M|tr 5 c and a by a factor 2. 



X 


dlogrd/dlogi (at 10 4 s) 


|M| tr 


2.5 xlO" 2 


c 


2.67 xl(T 3 


a 


-4.37 xlO -2 



cous outward transport of angular momentum that tends 
to make the disc larger. These effects are parameterized 
by Mtr, c and a respectively. In order to determine their 
influence on the disc radius variations, we have changed 
their values by a factor 2 (linear variation on a time scale 
of 10 4 s) for a steady disc in a hot state (Fig ||). The disc 
radius in the new steady state is unchanged, except in 
the case where c is modified; the short term influence of 
these variations are summarized in table || which shows 
the value of d log r/d log x, 10 4 s after the beginning of 
the variation, x being one of the three parameters. As can 
be seen, d log r jd log c is 5 to 10 times smaller than the two 
other quantities. This explains why ci/co has to be larger 
than 10 in order to obtain a superout burst, and why this 
ratio needs to be larger than 20 in order for the radius 
to decrease fast enough during a superoutburst and reach 
?"crito before the cooling wave starts. 

In principle, the disc shrinkage could be slower, or pos- 
sibly delayed if the mass transfer rate from the secondary 
were reduced, possibly as a result of a decrease of the il- 
lumination of the secondary. However, when the cooling 
wave starts, a decreases from «hot to ai ow in the outer disc 
regions, i.e. by a factor 5. c has to be kept fixed, since one 
requires the disc to remain eccentric. Because of the large 
change in viscosity when a cooling front starts to propa- 
gate, M tr would then have to be reduced by an extremely 



large factor; we made several numerical experiments with 
various prescriptions for a change in M tr (amplitude and 
possible delay with respect to the formation of a cooling 
front), and we were unable to find a case in which the 
minimum disc radius was not obtained shortly after the 
end of a superoutburst. 

4. Possible solutions 

4-.1. ER UMa supercycles 

The standard TTI model cannot reproduce the very short 
supercycles observed in ultra-low mass-ratio systems. Os- 



aki (1995a) reproduced RZ LMFs supercycle by ending 
prematurely the superoutburst assuming r cr i t o = 0.42a 
instead of the usual value of 0.35a (note that in the TTI 
framework, one assumes r cr i t o = 0.35a when modelling 
WZ Sge, even though its mass ratio is presumably smaller 
than in most SU UMa stars). In addition, the presence 
of superhumps after the end of the superoutburst implies 
that the eccentricity stops much later than the end of a 
superoutburst. The only ways for the disc to remain ec- 
centric after the end of a superoutburst are then: (i) tf is 
very long, or (ii) the accretion disc is always eccentric in 
ER UMa systems. 

In the first case, the transition time tf must be much 
larger than the one used here. However, as shown in Fig. 
||, the interval between a superoutburst and the next nor- 
mal outburst is far too long compared to observations. 
The second solution implies that the tidal torque is no 
longer responsible for superoutbursts, since c(r) = c\ is a 
constant. However, the presence of large outbursts is still 
possible, as the (modified) DIM predicts the alternation of 
narrow and wide outbursts (see paper I) for large enough 
mass transfer rates. If the mass transfer rate is constant, 
the light curve consists of one or two small outbursts sur- 
rounded by large ones, but small variations of the mass 
transfer rate can easily lead to ER UMa type light curves, 
provided that these variations show some regularity. Irra- 
diation of the disc and of the secondary can also account 
for the presence of long and short outbursts; Hameury et 
al. ( pOOO ) have included these effects in the standard disc 
instability model and produced light curves typical of sys- 
tems such as RZ LMi. 

We therefore conclude that ER UMa stars should be 
dwarf novae with a permanently eccentric accretion disc, 
thereby accounting for superhumps, and where the illumi- 
nation of the disc and the secondary star plays an impor- 
tant role. We thus predict that superhumps should exist 
at all phases of the supercycle of ER UMa stars; this ap- 
parently agrees with observations (Gao et al. 1999] ). 



4-- 2. Echo outbursts in WZ Sge stars 

WZ Sge stars have very long supercycles and superout- 
bursts. In addition, no normal outbursts between two con- 
secutive superoutbursts are observed. Some WZ Sge stars 
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also show echo outbursts at the end of the superoutburst: 
several small outbursts, spaced every tens of days, during 
which superhumps are still present. 



Smak ( 1993 ) deduces the mass transfer rate of WZ Sge 
from the luminosity of the hot spot. It is observed to in- 
crease at least by a factor 10 during superoutburst, and it 
decreases afterwards, remaining larger than the quiescent 
value during several tens of days. Such mass transfer rate 
fluctuation could result from irradiation of the secondary 
star by the white dwarf. Hameury et al. ( 200C| ) have shown 
that irradiation (including disc irradiation) could indeed 
account for some peculiarities of WZ Sge stars: for ex- 
ample, they reproduced the echo outbursts phenomenon 
without including the tidal instability in the DIM; how- 
ever, they did not reproduce a full cycle with long recur- 
rence times, and the echo outbursts they obtained were 
slightly too dim. 

The long recurrence times can be due to low alpha 
value (Smak 1993| ; Osaki 1995b), possibly due to a decay of 
the MHD turbulence that would lead to a time-dependant 
a— prescription (Osaki et al. 2001), thereby explaining the 
echo outbursts. A low viscosity could result from the sec- 
ondary being a brown dwarf (Meyer & Meyer- Hofmcistcr 



1999), but one would have to explain why the viscosity 



is so much lower in these systems as compared to other 
SU UMa systems which have comparable or even shorter 
orbital periods. Another possibility (Lasota et al 



1995 



1999; Hameury et al. 1997) is that WZ Sge stars are in 



a stable low state between superoutbursts, thus explain- 
ing the absence of normal outbursts. This requires a hole 
in the central regions of the disc, as a result of either a 
moderate magnetic field, or of evaporation. The superout- 
bursts would then be triggered by an externally imposed 
increase of the mass transfer rate, the long duration of the 
outburst and the large mass accreted onto the white dwarf 
being due to the irradiation of the secondary. They did not 
include the tidal instability and, as in the case of ER UMa 
stars, the presence of superhumps and late superhumps is 
explained if the accretion disc is always eccentric. If one 
combines these results with those of Hameury et al. ( [2000 ) 
on the echo outburst, one should be able to reproduce a 
WZ Sge star light curve with echo outbursts, which would 
however be slightly different than the observed ones. 

We have used our TTI model including irradiation and 
the presence of an inner hole for the accret ion di sc as de- 
termined by equation (7) of Hameury et al. ( 2000| ). We use 
the parameters of EG Cnc. The mass transfer rate from 
the secondary is assumed to be affected by irradiation ac- 
cording to: 

M tr = max(M 2 , 7 < M acc >) (4) 
where M% is the mass transfer rate that would be obtained 
in the absence of illumination, 7 is a constant in the range 
[0-1], and < Aiacc > is some average of the mass accretion 
onto the white dwarf, defined as: 



to being the actual time. This prescription differs slightly 
from that of Hameury et al. (1997), in that the average 
of the mass accretion rate is used instead of its value at 
time t$] this accounts for the fact that a fraction of the 
accretion luminosity is released at later times when the 
white dwarf cools. In the following, we have taken At = 
14 days. We start with a stable system in the low state 
with a mass transfer rate of 1.5 10 15 g.s -1 , in agreement 
with Smak ( |1993fl estimate for WZ Sge (2. 10 15 g s" 1 ). 
At t = 10 days, the mass transfer rate M2 is increased 
to 3 10 16 g.s _1 ; 20 days later, M2 is reduced to its initial 
value. Our results are insensitive on how M2 returns to its 
normal value, because illumination dominates in Eq. (^). 

This model predicts that the luminosity should in- 
crease just before the triggering of a superoutburst. The 
increase need not be large - one simply requires the mass 
transfer rate to exceed the critical rate for stability; this 
does t herefo re not contradict the observations of Ishioka 
et al. ( [2002 ) that no strong orbital hump due to a high 
mass transfer rate was detected during the early phase of 
WZ Sge 2001 outburst. 

We obtain a light curve with a very long superoutburst 
followed by 3 normal outbursts that last ~ 3 days and 
occur within two months of the main outburst (Fig. [5J). 
The disc radius remains always larger than 0.35 a, and 
is always in an eccentric state; permanent superhumps 
are thus expected. However, the interval between normal 
outbursts is twice longer than for EG Cnc echo outbursts, 
which are brighter than observed. We did not obtain a 
better fit with ER UMa light curves; in view of the very 
crude assumptions made to derive the mass transfer rate 
from the secondary under the influence of illumination, 
this is not really a surprise. 

5. Conclusion 



Hellier (2001) suggested that several limitations of the 
TTI model could be overcome if, in ultra-low mass ra- 
tio systems, the thermal instability ends while the disc is 
still eccentric. This would explain the late superhumps, 



and Hellier (2001) also suggested that this would account 



< M a 



>= 



-(t -t)/At dt 



(5) 



for the very short ER UMa supercycles and the WZ Sge 
stars echo outbursts. 

We do find that the decoupling of the two instabilities 
can happen much more easily in ULMR systems than in 
other SU UMa stars. However, the tidal instability stops 
almost immediately after the end of the thermal-viscous 
instability (delay less than a day) , because the disc shrinks 
very rapidly after the departure of a cooling wave. We esti- 
mate that it takes only a few days for the disc to return to 
a circular shape after the cessation of the tidal instability. 
This is still insufficient to explain late superhumps that 
are found tens of days after the superoutburst in some 
systems. 

We propose instead that the discs in ER UMa systems 
are permanently eccentric, and that they should therefore 
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Fig. 5. Response of a steady accretion disc to a mass 
transfer rate variation by a factor 20. We use the parame- 
ters of EG Cnc, include irradiation, the presence of a hole 
in the disc, and the existence of a tidal instability. Initially 
Mi = 1.5 10 15 g.s -1 ; at * = 10 M 2 increases up to 3.0 10 16 
g.s -1 , and at t — 30 returns to its quiescent value. The top 
panel shows the mass accretion rate onto the white dwarf 
(solid line) and the mass transfer rate from the secondary 
(dashed line); the other panels show the outer disc radius, 
its mass and visual magnitude 

have permanent superhumps; there are observational in- 
dications of this (Gao et al. 1999Q . The short supercycle 
is then due to either an alternation of narrow and wide 



outbursts as in SS Cyg (see Buat-Menard et al. 2001 ), or 
to the effects of irradiation (see Hameury et al. 20001) . Fur- 
ther observations of late superhumps should confirm their 
permanence in ER UMa stars. 



As in Hameury et al. (1997), we propose that WZ Sge 



transfer rate, and echo outbursts are triggered while the 
mass transfer rate is still high after the superout burst. 
However, simulated light curves are still somewhat dif- 
ferent from observations; this might be due to our crude 
treatment of the secondary illumination. 
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